Dynamics and chaos in the unified scalar field Cosmology 
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We study the dynamics of the closed scalar field FRW cosmological models in the framework of 
the so called Unified Dark Matter (UDM) scenario. Performing a theoretical as well as a numerical 
analysis we find that there is a strong indication of chaos in agreement with previous studies. We 
find that a positive value of the spatial curvature is essential for the appearance of chaoticity, though 
the Lyapunov number seems to be independent of the curvature value. Models that are close to 
flat (fc — > + ) exhibit a chaotic behavior after a long time while pure flat models do not exhibit 
any chaos. Moreover, we find that some of the semiflat models in the UDM scenario exhibit similar 
dynamical behavior with the A cosmology despite their chaoticity. 

Finally, we compare the measured evolution of the Hubble parameter derived from the differential 
ages of passively evolving galaxies with that expected in the semiflat unified scalar field cosmology. 
Based on a specific set of initial conditions we find that the UDM scalar field model matches well 
the observational data. 
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1. INTRODUCTION 

There is by now convincing evidence that the avail- 
able high quality cosmological data (Type la supernovae, 
CMB, etc.) are well fitted by an emerging cosmologi- 
cal model, which contains cold dark matter to explain 
clustering and an extra component with negative pres- 
sure, the vacuum energy (or in a more general setting the 
"dark energy"), to explain the observed accelerated cos- 
mic expansion (Refs. [3, [|J and references therein). Due 
to the absence of a physically well- motivated fundamen- 
tal theory, there have been many theoretical speculations 
regarding the nature of the above exotic dark energy. In- 
deed, many authors claim that a real scalar field which 
rolls down the potential U(cj>) 0, 0, 1, S 0, S, ©, El could 
resemble the dark energy. 

On the other hand, the global dynamics, in models 
with dark energy, become more complicated than in mat- 
ter dominated models. Within this framework, a serious 
issue is the existence (or nonexistence) of chaos in the 
scalar field cosmology. This is very important because 
chaotic fields can provide possible solutions to the cos- 
mological coincidence problem (the question why Q m and 
SI A are of the same order at the present time) as well as 
to the problem of uniqueness of vacua. The first paper 
on chaos in a closed Friedman Robertson Walker geom- 
etry coupled to a scalar field (FRW with a scalar field 
potential U(<f>) cx (f> 2 ) was written by Page [ll[. This au- 
thor, suggested that the dynamical behavior of this non- 
integrable model is similar to the chaos appearing in er- 
godic systems. Since then, the problem of chaos in FRW 
scalar field cosmologies has been investigated thorou ghly 
in several papers [lfl 0, El M El El OS EE M SlMf • 
The dynamical properties of such cosmological models 
are met in dynamical systems exhibiting chaotic scatter- 
ing. However, in this kind of studies there is a crucial 



question which has to be addressed: how can we esti- 
mate the amount of chaos? It has been shown that in 
chaotic scattering problems classical tools, like the Lya- 
punov characteristic number, are unable to provide a 
clear answer whether an orbit is chaotic or not (22J. 

The aim of this work is to investigate the dynamical 
properties of a family of cosmological models based on a 
scalar field potential which appears to act both as a dark 
matter and dark energy [1 [23| . In order to measure the 
amount of chaos we follow the notation of Ref . [22j that 
was able to give useful answers to their scattering prob- 
lem. These authors used an index of chaoticity which 
is capable to detect stochastic areas in scattering prob- 
lems. To this end, we attempt to compare our theoretical 
predictions with observational data by utilizing the mea- 
sured evolution of the Hubble parameter [24| • 

The structure of the paper is as follows. The basic the- 
oretical elements of the problem are presented in section 
2. In particular, in the first subsection we study the dy- 
namical system analytically and in the second subsection 
we consider various chaoticity indices. Section 3 outlines 
the numerical results of the present study together with a 
thorough discussion, and finally we draw our conclusions 
in section 4. 



2. THEORETICAL ELEMENTS 



2.1. Robertson- Walker cosmology coupled with a 
scalar field. 



Observationally, the universe at large scales appears 
to be in general homogeneous and isotropic. Such a uni- 
verse can be described by the Robertson- Walker (RW) 
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line element: 

, 2 , 9 n, , d 2 r + r 2 d 2 6 + r 2 sin 2 6d 2 cf> 

ds 2 = -dt 2 + a 2 (t) — 2 (1) 

v i (1 + ^-) 2 

where a(i) is the scale factor of the universe and the 
sign of k determines the geometry of the space. The 
RW cosmology has been supplemented in recent years by 
the theories of inflation, dark matter and dark energy. 
Many of these complementary theories are based on the 
presence of one or more scalar fields. In this work we 
have ignored both the coupling of the scalar field </> to 
other fields and quantum-mechanical effects. Doing so, 
it is straightforward to associate the gravitational RW 
universe with the scalar field by utilizing the so called 
stress-energy tensor (see Ref. [25[): 
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I • 



T 
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J 2 
TH = 



! + U{4>) 

U{cf>) 



for i = j 
for fjt^u 



(2) 



where the overdot denotes derivatives with respect to 
time, U(<fi) is the potential of the scalar field, the Latin 
indices indicate the space coordinates, while the Greek 
indices the whole spacetime. 

Using now the Einstein's field equations [33j 



(3) 



and the equation of state through conservation of the 
stress-energy: 



T — 



the corresponding Friedmann equations become: 
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(5) 
(6) 



while the equation of motion for the scalar field takes the 
form: 



. + 3^+^ = 
a a<p 



(7) 



The above set of differential equations describes the 
evolution of the Robertson- Walker cosmology when grav- 
ity is coupled with a scalar field. It can be proved (see 
Ref. [ll], [26j]) that these equations of motion can also 
emerge from the Lagrangian: 



L 



-3aa + Ska + a 



(8) 



A serious problem that hampers the straightforward use 
of such an approach is our ignorance of the form of the 
potential U {(f)). Note that in the literature, because of the 
absence of a physically well-motivated fundamental the- 
ory, there are plenty of such potentials (for a review see 
|2g |) which approach differently the nature of the scalar 
field. For example, if the main term in the potential 
U((j)) is 4> n then the energy density of the scalar field is 



oc a 



-6n/(n+2) 



Therefore, it becomes evident that for 



n = 2 or n = 4, the corresponding energy density behaves 
either like non relativistic or relativistic matter [25| . 

In a recent work [23j, the authors used the so called 
Unified Dark Matter (UDM) scenario in order to param- 
eterize the functional form of the potential 



U{4>) = a cosh 2 (c 4>) + 



(-2 



(9) 



From a cosmological point of view, this "cosmic" fluid be- 
haves both as a dark energy and dark matter. However, 
at an early enough epoch the fluid evolves like radiation 
[27l ]. Also it has been found [23[ that in a spatially flat 
RW model (k = 0) the UDM fluid predicts the same 
global dynamics as the A cosmology does. 

The real constants c\ and C2 in Eq. ([9]) obey certain 
restrictions because the square of the mass of the scalar 
dJIW > for = and furthermore, 
U(0) > 0, otherwise the scalar 



field is equal to 



~5S£ 

according to Ref. [2S 
field mass could get negative values. Explicitly these re- 
strictions mean for c\ and C2 that: 



ci > 
ci > -c 2 



(10) 
(11) 



Changing variables from (a, < 
tions 



to (x, y) using the rela- 



x = A a 3/2 smb. (c 0) 

y = A o?l 2 cosh (c (j>) (12) 

with A 2 = 1/c 2 = 8/3 the Lagrangian ([S]) is written: 



L 



{x 2 + -c 2 x 2 ) 



(y 2 + -(ci+c 2 ) y 2 ) 



^k{ 



y 



(13) 



Hence, it is obvious that in the new coordinate system our 
problem is described by two coupled oscillators. Notice 
that the oscillator on the y axis is hyperbolic due to the 
restriction (fTTj) . but the oscillator on the x axis is either 
hyperbolic, if c 2 > 0, or elliptic, if c 2 < 0. Also note that 
x and y are constrained due to Eq. (|12p which gives: 



y > 0, y > \x\ 
because the scale factor 

-2\ 1/3 



y 



A 2 



(14) 



(15) 
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is physically meaningful only if it is greater or equal to 
0. 

As can be seen from the above analysis there are two 
cases of interest: 

1. c 2 > 

2. c 2 < 

and we investigate them in the subsequent subsections. 
Finally, it is obvious that for a flat cosmological model 
(k = 0) the dynamical system is fully integrable. 



2.1.1. Case c 2 > 

In this CcLSG, cLS mentioned before, both oscillators are 
hyperbolic. From the Lagrangian ()13|) we find the Hamil- 
tonian: 



n= l - [(pi - tujx 2 ) - ( Py - 



1 r 



34/3^2 _ ^2)1/3 



(16) 



where p x = x, p y — —y denote the canonical momenta 
and uj\ = |c2, = |( c i + c ri) are the oscillators' fre- 
quencies. 

The Hamiltonian (I16|) gives the following equations of 
motion: 



Py 



dH 

dp x 

on 

dp y 
dH 
dx 
dH 
dy 



Px 



-Py 



-CJ? + 



(y 2 - x 2 ) 2 / 3 
3^ 3 k 



x (17) 



(y 2 — x 2 ) 2 / 3 



II 



From Eqs. (jTTJ) and the restrictions (TT3]) we find that the 
equilibrium point of the system is: 



3 1 / 4 k 3 / 4 
3/2 ' 



Py = Px = X = 



(18) 



when k > 0. If k < 0, then there are no equilibrium 
points which obey the restrictions (|14j) . 



2.1.2. Casec 2 <0 

In this case we get one hyperbolic and one elliptic os- 
cillator and the Hamiltonian becomes: 



* = £[(pS + a**V0£-^V)] 



3 4 / 3 % 2 - x 2 ) 1 ' 3 



(19) 



where io\ = — \c%. 

The Hamiltonian (JTHJ) now provides a somewhat differ- 
ent set of equations than before: 



Px 



Py 



m 

dp x 
&H 

Spy 

an 

dx 

m 

dy 



Px 



-Py 



-w + 



3V3fc 



(>J 2 



2-2)2/3 



(20) 



3V3jfe 



(y 2 - x 2 ) 2 ' 3 



From Eqs. P0|) . and by taking into account the restric- 
tion (fT4]) , we find that the equilibrium point of the system 



is again: 



3i/4p/4 
3/2 ' 



Py = Px = X = 



(21) 



when k > 0, which is the same as in the previous para- 
graph (see Eq. (jX8j) ) . If k < then there are again no 
equilibrium points. The physical interpretation of the 
non existence of equilibrium points when the universe is 
flat or has negative curvature is that the universe can- 
not collapse therefore there is no other flow in the phase 
space than the one which starts with a = and tends to 
a — > 00. 

For both cases of C2 at the equilibrium point the eigen- 
values and the corresponding eigenvectors are: 



±a/-|ci (Px,Py,x,y) = (±\Pfi, 0,1,0) (22) 
±Vci + c 2 (p x ,Py,x,y) = (0, ±Vc~i + c 2 ,0, 1)(23) 

Due to the restriction (|10|) the eigenvalues of (|22|) are 
imaginary, but due to the restriction (fl"Tj) the eigenvalues 
of (|23p are real, so we have a saddle-center instability. 
The saddle shape is displayed on the y, p y plane. No- 
tice that, although the existence and the position of the 
equilibrium point depend from the curvature's value, the 
eigenvalues and the corresponding eigenvectors for the 
equilibrium point are independent from the curvature. 

2.2. Chaotic indicators and scattering 

The detection of chaos (see [28] for a review) is a topic 
for which many tools have been developed. These meth- 
ods can be separated in two main categories. Tools which 
exploit the properties of the deviation vector and tools 
that are based upon a frequency analysis. In this paper 
we use tools of the first category. The deviation vector £ 
is defined as the solution of the variational equations of 
the system. 

A classical tool for the detection of chaos is the Lya- 
punov Characteristic Number (LCN): 



LCN = lim X (t), 

t^OO 



(24) 
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where: 



x[) t c(o) 



(25) 



is the "finite time LCN" . The orbits with LCN > are 
called chaotic, while the orbits with LCN=0 are called or- 
dered. A similar useful quantity is the stretching number 
|29j |. also called the Lyapunov indicator (30j |: 



sn(t) = In 



£(t + dt) 



(26) 



In scattering problems when orbits escape from the re- 
gion where they behave stochastically, and then get far 
away from this region, they diverge from nearby orbits 
linearly. Therefore the LCN is zero. The latter behav- 
ior is also a characteristic of regular orbits. Hence if we 
entrust entirely the LCN in scattering problems we may 
characterize as regular those orbits, that in a certain re- 
gion exhibit stochastic behavior. However in our paper, 
in order to study chaos we use also the evolution of the 
stretching number, as done in Ref. This method is 

based on the idea that when an orbit in a chaotic scat- 
tering problem passes through a stochastic region of the 
phase space the average value of stretching varies around 
a positive number, even if afterwards it tends to a value 
near zero. 



3. NUMERICAL RESULTS. COMPARISON 
WITH OBSERVATIONAL DATA 

In the present section we investigate numerically the 
system of the RW coupled with a scalar field when the 
energy 7i is zero. Our investigation refers mainly to the 
case C2 > following the notations of [23|. Note that in 
the appendix we present results for the case C2 < 0. If 
we introduce the noncanonical transformation (k > 0) 



P 



x = 4 X y = fc 3 / 4 r 



(27) 



then the Hamiltonian (|16p gives 

k 3/2 

n = — — x 



(P| - cofX 2 ) - (P 2 _ ^y 2 ) - 3 4 / 3 (F 2 - X 2 )^ 3 ](28) 



The Hamiltonian (TT9|) for C2 < transforms in the same 
way. The only difference is that we have to change 
the sign of ojf in (]28[) . So as regards the parameter k 
for TL = we may consider only one value of k (here 
k = 1CP 3 , hereafter called semi-fiat case), because the 
transformation (|2T|) shows that any case of k can be 
transformed to another. However, we must stress that 
the time rescales as k 3 ^ 2 by the above transformation. 
For other values of the curvature k the value of t must be 
multiplied by k 3 ! 2 . The reason we choose the semi-flat 



value of k is that the analysis of the recent observations 
of the Cosmic Microwave Backround (CMB) anisotropics 
indicates that the spatial curvature of the universe takes 
a quite small positive value (fifc w 0.02) very close to 
zero [H • Our aim in this work is to investigate whether a 
small variation of the curvature around its nominal value 
(k = 0) could affect the global dynamics of the universe. 
We do not present the cases where k < because of the 
recent observational data. Furthermore, in the case of 
the flat universe the two oscillators are uncoupled so the 
system is integrable. 

In the case of two coupled hyperbolic oscillators we 
select c\ = 1 and ca = 1 for our parameters, which means 

that the potential © is simply £/(</>) = cosh 2 (y§0) + L 
Even such a simple parameterization produces a variety 
of orbits. In order to show this variety when k = 10~ 3 we 
choose five representative orbits. The initial conditions 
of these orbits are: 1) y = 0.0001, x = 0, p x = 0.006, 2) 
y = 0.0001, x = 0, p x = 0.03, 3) y = 0.0001, x = 0, p x = 
0.04, 4) y = 0.0001, x = 0, p x = 0.041, 5) y = 0.0001, 
x = 0, p x — 0.099. For all these orbits p y is evaluated by 
using Eq. (fl"6| with H = and by choosing the negative 
sign in front of the square root of p 2 . 

Figure [IJi shows the scalar field as a function of the 
logarithm of the scale factor log 10 a for our five orbits 
and in the Fig. [TJ) the logarithm of the scale factor as a 
function of the logarithm of time log 10 t. We see that the 
orbits 1, 2, 3 refer to universes that recollapse, therefore 
these orbits describe examples of closed universes. On 
the other hand the orbits 4 and 5 refer to universes that 
expand to infinity, which means that even with a positive 
curvature we get cosmological models which expand for 
ever (see Ref. 



31 



It is obvious that as the curvature k 
increases the orbits which describe closed universes will 
reach bigger values of the scale factor before they col- 
lapse, which also suggests that they will collapse after a 
longer time. It is interesting to mention that for those 
orbits which are in between 4 and 5 the scale factor of 
the universe initially goes as oc i 1 / 2 , then for a certain 
interval of time as oc t 2 ^ 3 , and finally accelerates expo- 
nentially. In order to scale the dynamic time t of the 
system to the real time t\j in Gyears, we fit our numeri- 
cal results for the orbit 5 to the observational data. That 
way we find the relations tjj ~ 6.1i Gyrs and ajj = 3a. 

Although, the (a, </)) coordinate system is the physi- 
cally straight forward for the present problem, the (x, y) 
system gives a more clear insight to the problem. Thus in 
Fig. Wfr w e present the projections of the orbits on the y, 
p y plane, while in Fig. (2b we give the projections of the 
orbits on the x, p x plane. In the y, p y plane the orbits 1, 
2, 3 start from y = and finally return to y — 0. This is 
expected simply because these orbits model recollapsing 
universes a-* 0. The orbits 4 and 5 describe expanding 
cases where y, like the scale factor, tends to infinity and 
p y — ► — oo. The form of our orbits in the x, p x plane 
are similar, except for the orbit 4 which, instead of giv- 
ing x — > co returns to x = with p x — > 0. This means 
that the orbit 4 is essentially driven by the y 7 p y oscil- 
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FIG. 1: (a) The scalar field cj> vs. the logarithm of the sc 
logarithm of time log 10 1 for 5 orbits in the semi-flat case (A; 



(a) 




FIG. 2: (a) The projections of the orbits on the y, p y plane, 
in the semi-flat case (k — 10~ 3 ). 

lator which gains all the energy while the x, p x motion 
ceases to exist. Thus, the case of orbit 4 shows explicitely 
that the two coupled oscillators exchange energy. The en- 
ergy exchange takes place due to the coupling term. The 
question which arises now is if this exchange is done in 
a regular way or if it is done stochastically. To address 
this question we used the variational equations of (fT7|) to 
calculate the evolution of the stretching numbers and of 
the finite time LCN "x" through time. 

Doing so, in Fig. [3] we present our results: in Fig[3h 
is shown the stretching number (hereafter "sn") versus 
the logarithm of time log 10 1 and in Fig|3h> the logarithm 
of x (l°gio x) versus the logarithm of time log 10 1. The 
dashed line in the FigJHjD represents the inclination — 1 



i factor log 10 q. (b) The logarithm of the scale factor vs. the 
lO" 3 ). 
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(b) The projections of the orbits on the x, p x plane for 5 orbits 



which is followed by the regular orbits as the LCN tends 
to 0. This dashed line helps us to separate the regu- 
lar from the chaotic orbits, because the regular orbits 
always evolve along a parallel line, though usually with 
some variations around the line, while the chaotic orbits 
change after a certain time point their inclination from -1 
to as the \ tends to a finite value. For all the considered 
orbits Fig [3^ the "sn" stays for most of the time above 
and this indicates that for all five orbits x should be 
positive, therefore there is a strong indication that these 
orbits are chaotic. This conclusion seems to agree with 
the results which we see also for the orbits 4 and 5 in Fig. 
[3Jd, but the orbits 1, 2 and 3 are not so clear cases because 
they recollapse. In particular, the expanding to infinity 



6 




O 



1 - 



- 




FIG. 3: (a) The stretching number "sn" versus the logarithm of time log 10 t. (b) The logarithm of the finite time LCN, 
logiox, versus the logarithm of time log 10 t for 5 orbits in the semi-flat case (k = 10 -3 ). The dashed line in (b) represents the 
inclination —1 which is followed by the regular orbits. The initial, normalized to unity, deviation vector for all the orbits is 
1, 1, 1) and the step of integration is dt = 10 -5 . 




FIG. 4: The Hubble parameter versus the scale factor of the 
universe. The points corresponds to the observational data. 
The gray dashed line represents the A cosmology. The solid 
black line corresponds to the the UDM model with k = 10~ 3 
and initial conditions x — 5.6 x 10 -7 , y — 0, p y = 0.09. Note 
that the scale factor scaled to unity at the present time. 



orbits 4 and 5 appear to be chaotic. Their \ ( see Fig. EJd) 
tends to a finite nonzero value « v2 and the "sn" is on 
the average positive. This finite \ value is the eigenvalue 
V c i + c 2 (|23p of the unstable manifold, which emanates 
from the equilibrium point, for c\ = c% = 1. Notice that 



this eigenvalue does not depend on the curvature, hence 
for any positive curvature in the case under study we 
expect that the chaotic orbits will have the same LCN. 
This expectation was numerically attested for several or- 
ders of magnitude of the curvature's value in the case 
C\ = C2 = 1. This is true if the orbits do not recollapse. 
If the orbits recollapse before the x indicator can decide 
whether they are chaotic, then the "sn" is the only indi- 
cator that can suggest the chaoticity of these orbits. In 
the case of the orbits 1, 2 and 3 (Fig. [3^) the "sn" indeed 
indicates that these orbits are chaotic. 

The recollapse for the orbits 1,2 and 3 in Figs. [5^,b 
can be observed by an abrupt change of "sn" and \ in- 
dicators respectively. At the recollapse both indicators' 
values start to fall, abruptly for orbits 1 and 2, marginally 
for orbit 3, before they steeply rise and become infinite. 
The fall of the indicators at the recollapse is explained 
by the fact that when an orbit falls to the anomaly a = 
the position variables tend to stabilize, thus the devia- 
tions of the position variables shrink and the deviation 
vector's norm gets smaller. As long as the value of the de- 
viation vector's norm decreases, the "sn" becomes more 
and more negative and therefore the \ value decreases. 
However, as the moving point gets nearer to the anomaly 
the momentum variables start to grow and their devia- 
tions expand. This means that the deviation vector's 
norm increases, thus the "sn" increases and the \ value 
rises. Because the momentum variables grow abruptly 
the rise of "sn" and \ is steep. 

We should state here that the anomalies produced by 
the recollapse could be avoided by a logarithmic trans- 
formation of the time, as done in the mixmaster problem 
(for a review see [28|). However, we do not discuss here 
the usefulness of such a transformation, because the re- 
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cent high quality cosmological data show that the uni- 
verse is expanding accelerated to infinity and does not 
recollapse. 

We emphasize that the final value of x does not depend 
on the curvature's value k. Hence cosmological models 
for any positive value of k reveal a chaotic behavior in 
contrast with the pure flat models (k = 0). As a special 
case, semiflat models (k — 10~ 3 ) with the characteristics 
close to orbits 4 and 5 predict an overall dynamics which 
is close to the A cosmology [13], despite the fact that 
for these UDM cosmological paradigms there is a strong 
indication of a chaotic behavior! Indeed, orbit 5 starts 
with a state which is similar to the radiation epoch (a oc 
i 1 / 2 ), then evolves to the matter epoch (a oc i 2 / 3 ) and 
finishes with the exponential growth of the universe 

On the other hand, it is interesting to mention that 
the interplay between the values of the initial conditions 
could yield semiflat UDM cosmological models which pre- 
dict well the observed evolution of the Hubble parame- 
ter. As an example, in Fig. 4 we plot the measured 
Hubble parameter together with a UDM scenario (solid 
line) which is formed by the following initial conditions: 
x = 5.6 x 10~ 7 , y = 0, p y = 0.09. It is obvious that 
this model reproduces, even better than the A cosmol- 
ogy (dashed line), the observational data. Note that the 
chaotic behavior of this UDM paradigm is the same with 
that found in the case of orbit 5. In a subsequent paper 
we shall investigate statistically the parameter space for 
models with k > which fit best the available observa- 
tional data. 

The above result is novel because within the frame- 
work of the scalar field cosmology we can find semiflat 
models with specific initial conditions which reveal chaos 
and at the same time provide an evolution of the Hub- 
ble parameter close to the observed one. In any case, we 
expect that the exact value of the curvature will be fi- 
nally measured by the next generation of the CMB data 
(PLANCK observatory) and therefore, we cannot pre- 
clude observational surprises regarding the curvature of 
the universe. 



4. CONCLUSIONS 

In this work we investigate the dynamics of the closed 
scalar field FRW cosmological models in the framework 
of the so called Unified Dark Matter scenario. We find 
that for the closed geometry there is a strong indication 
of chaos in agreement with previous studies and that the 
LCN of these models is independent from the curvature. 
In particular, we find that there are semiflat cosmological 
models with specific initial conditions in which there is 
a clear indication for a chaotic behavior and at the same 
time the corresponding dynamics is almost the same as 
in the A cosmology. We verify this by combaning the 
measured evolution of the Hubble parameter with that 
expected in the scalar field cosmology (with specific ini- 
tial conditions) and we find a very good agreement. If 
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FIG. 5: (a) The Poincaree section x — with p x > 0. (b) 
The orbits' projections on the x, p x plane for 4 orbits, (c) 
The development of log 10 a as a function of the logarithm of 
time log 10 t with k = 10 -3 . 
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that is the case, then the chaotic fields may provide an 
alternative theory for the solution of the cosmological co- 
incidence problem. 
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APPENDIX A: CASE c 2 < 

In this kind of dynamical studies, it is well known (for 
details see [28| ) that one of the best procedures to find 
chaos is based on the so called Poincare sections. Unfor- 
tunately, in the case of C2 > the Poincare sections are 
not applicable, due to the fact that the orbits do not ex- 
hibit recurrences. On the other hand, for C2 < Poincare 
sections can be obtained because of the elliptic oscillator 
on x, p x plane and that only if the frequency of the el- 
liptic oscillator is much higher than the frequency of the 
hyperbolic one. Hence in order to achieve the latter we 
use here c\ — 10 4 + 1 and = 1 — 10 4 . We mention 
that the case ci < is mainly of mathematical interest, 
because, as shown in Ref. [231 ] , C2 is meaningful only if it 
is not negative. 

In Fig. are shown Poincare sections for k = 10 -3 . 
The surface of section is for x = and p x > 0. From these 
Poincare sections we conclude that the present dynamical 
system's phase space is separated into 4 subspaces with 
specific directions of flow. Each subspace has a typical 
type of orbits, therefore we need only 4 appropriate initial 
conditions in order to represent the whole section's phase 
space. The initial conditions we choose are 1) y = 0.001, 
p y = -0.0072, x = 0, 2) y = 0.001, p v = -0.0068, 
x = 0, 3) y = 0.011, p y = 0.0072, x = 0, 4) y = 0.011, 
p y = 0.0076, x = 0. Note that in this subsection p x is 
calculated from Eq. (|19|) with H = and by choosing 
the positive sign in front of the square root of p\ . The 
direction of the flow for orbits of type: 

1. is along the y direction. 

2. is along the p y direction. 

3. is along the negative y direction. 

4. is along the negative p y direction. 



The equilibrium point in the Poincare section is at the 
center of the cross formed by the representative orbits, 
as given by Eqs. I|21[) . The elliptic oscillator can be seen 
in Fig. [SJd where we give the projection of each orbit 
on the x, p x plane. As expected, the orbits on the x, 
p x projection move almost on ellipses whose axes have 
a ratio ^ ~ -j^. The numbers in Fig. [SJs written on 
the x = axis indicates how the four orbits are ordered 
according to their p x axis length. We mention that the 
length of each axis varies as a function of time. 

From the dynamical point of view, which we get from 
Figs. [S^i, b, we now move on to the more physical point 
of view of Fig. [5};. In this figure we give the development 
of log 10 a as a function of the logarithm of time log 10 t. 
In these figures we see orbits of the following variety. 
Models: 

1. describe a cosmological model which starts from 
a — > and ends with a — > oo after an inflection 
point. 

2. describe a cosmological model which starts from 
a — > and collapses. 

3. describe a cosmological model which starts from a 
finite a and ends with a — > oo. 

4. describe a cosmological model which starts from a 
finite a and collapses after an inflection point. This 
infection point represent the passing of the scale 
factor from a > to a < 

Models of type 1 are the only ones that could be com- 
patible with the standing cosmology, but the scale factor 
at first increases very slowly with t and in the end it in- 
creases superexponentially. As for the chaoticity of the 
orbits in the C2 < case, the indexes \ an d "sn" indicate 
that these orbits are chaotic. 

We studied also some cases with H ^ for both 
C2 > and C2 < 0. The whole phase space in all cases 
is dependent on the position of the equilibrium point, 
because the asymptotic manifolds emanating from the 
equilibrium point determines the flow in the phase space 
and therefore determines which models will collapse and 
which shall expand to infinity. As can be seen from Eqs. 
([T5)) . (f2"T) the equilibrium point is independent from H, 
so the main characteristics of the 7i ^ cases are similar 
to those with Ti — 0. 
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